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PLh ■ Abstract 

The paper considers instantly coalescing, or instantly annihilating, systems of 
one-dimensional Brownian particles on the real line. Under maximal entrance laws, 
the distribution of the particles at a fixed time is shown to be Pfaffian point processes 
closely related to the Pfaffian point process describing one dimensional distribution 
of real eigenvalues in the real Ginibre ensemble of random matrices. As an applica- 
tion, an exact large time asymptotic for the n-point density function for coalescing 
• particles is derived. 

in 

q ■ 1 Introduction and summary of main results 

O ' 

The study of single species reaction diffusion systems A + A — > A (coalescence) and 
A + A — > (annihilation) originated in non-equilibrium statistical mechanics (see |13j). 
but has now a large mathematical literature (see, for example, [I], [6], [7], [I]). In 



one dimension the systems exhibit strongly non-mean field behaviour due to correlation 
effects. In this paper we give several examples showing that this correlation structure can 
be encoded algebraically in a Pfaffian structure. Note that the embedding of annihilating 
random walks as domain boundaries for a Glauber model makes Pfaffian formulae quite 
reasonable due to the free fermion structure of the Glauber model (see Felderhof [8]). 

We examine the asymptotics for the n-particle density function for (instantly) 
coalescing Brownian motions on R defined by 

P [there exist particles in dxi, . . . , dx n at time t] = p[ n \x\, . . . , x n ) dx\ . . . dx n . 
In [12] we showed, for n > 1, to,L > and for a variety of initial conditions, the bounds 

< Cl (n,L,t ) < - n( „C (ai '" — < c 2 (n,L,t ) < oo, (1) 

t 2 4 Ul<i<j<n\ x i ~ X J\ 
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for t > to and |acj| < Lt 1 / 2 , where the constant c\ will depend also on the initial condition. 
The non-linear factor in the power of t illustrates the non mean-field behaviour due to 
correlations. 



In this paper we show that, under the maximal entrance law, the true asymptotic holds 
in dU as t — > oo, and identify the limiting constant as the Pfaffian of a certain matrix. The 
maximal entrance law corresponds intuitively to starting with every point occupied, and 
can be constructed as the limit of initial Poisson distributions with increasing intensities. 
This initial condition is natural since, as explained in section 12.31 started from a large 
class of other initial conditions the distributions at time t become close, as t — > oo, to 
those of the maximal entrance law. 



Theorem 1 Under the maximal entrance law for coalescing Brownian motions, 



sup 

\x t \«t^ 2 



p? n \ xl ,...,x 2n ) 



n(2n-l) 

t- n -^-^ n 



(47r)- n/2 P/(J (2n) (0)) 



l<i<j<2n\ X i X j\ 



as t — ^ OO, 



where \x{\ « t 1 ^ 2 means that we may take the supremum over any positions (xi(t)) 
provided that supj \xi(t)\t~ 1 / 2 as t — > oo, and where Pf(j( 2n \<f))) is the Pfaffian of 
the In x 2n anti- symmetric matrix J^ 2n \4>) with entries 



4f ) w = (-ir 1 



(i - l)!(j - 1)! dx i+ i 



d i+j ~ 2 <±> 

1(0) forl<i< j < 2n, 



where <p(z) = zexp(— z 2 /A). Under the analogous maximal entrance law for annihilating 
Brownian motions, the same limit holds with (47r) -ra / 2 replaced by (64-7r)~ n / 2 . 



The presence of the Pfaffian in this asymptotic is a reflection that under the maximal 
entrance law the particle positions, at a fixed time, form a Pfaffian point process (see 
the start of section [3] for a definition) . 



Theorem 2 Under the maximal entrance law for coalescing Brownian motions, the par- 
ticle positions at time t form a Pfaffian point process with kernel t~ l l 2 K(xt~ l l 2 ,yt~ l l 2 ), 
where 

-F"(y - x) -F'(y - x) 

F'(y-x) sgn(y - x)F(\y - x\) 



K(x,y) 



and F(x) = 7r _1//2 e~ z2 ^dz. (Here sgn(z) = 1 for z > 0, sgn(z) = —1 for z < 
and sgn(0) = 0.) Under the analogous maximal entrance law for annihilating Brown- 
ian motions, the particle positions at time t form a Pfaffian point process with kernel 
\t- x l 2 K{xt- x l 2 ,yt- x l 2 ). ' 



The annihilating versions of Theorems Q] and [2j that is for (instantly) annihilating parti- 
cles, can be deduced from the thinning relation that connects coalescing and annihilating 
systems (see section [27T]) . 
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Many probabilities for the fixed t distributions are given by formulae using Pfaffians, 
and there are many places to start when proving these formulae. We choose to start by 
considering the following basic fact for product moments for annihilating systems, from 
which we will deduce all the other Pfaffian fromulae. 

Theorem 3 Consider the product moments for annihilating Brownian motions, defined 
by 

for bounded measurable g, where (x\, . . . , x n ) lists the initial positions of the annihilating 
Brownian motions on R, and [X\ : i G It) list the positions of any particles that remain 
at time t (and an empty product is taken to have value 1 ). Then for x\ < X2 < . . . < X2 n , 
the even moments m[ 2n \x) are given by 

mf n \ Xl , . . . , x 2n ) = Pf (mf\x h Xj ) : 1 < i < j < 2n) (2) 

where the right hand side is the Pfaffian of the 2n x 2n anti- symmetric matrix with 

(2) 

entries ml (xi,Xj) above the diagonal. 

Note these Pfaffians are in variables that determine the initial conditions, allowing us 
to use p.d.e. methods to characterize these moments. Indeed, the product moments 
satisfy a closed system of heat equations (with suitable boundary conditions), and we 
will verify Theorem [3] by simply checking that the Pfaffian uniquely satisfies this system. 
Markov time-reversal duality (see section 12. 2\i then immediately implies that certain 
empty interval formulae 

P [the intervals (01,02), (03, ct4),...,(a2m-lj «2m) are empty at time t] 

for coalescing systems are given by a Pfaffian, where the Pfaffian is now in the variables 
«i < 0-2 < ■ ■ ■ < a2m that determine the end points of the target intervals. This quickly 
leads to the identification of the Pfaffian point process kernel K(x,y). 

We concentrate on Brownian particles but, as we note later, we expect many of our 
Pfaffian formulae to hold for a large variety of spatial motion processes, and the Pfaffian 
structure seems to arise from two basic underlying mechanisms: linearly ordered particle 
motion and instantaneous reactions. 



1.1 Relation between annihilating Brownian motions and the real Gini- 
bre ensemble of random matrices. 

The Pfaffian point process defined in Theorem [2] has been originally discovered in the 
context of random matrices Namely, consider real Ginibre ensemble [10] defined by the 

1 We are grateful to the anonymous referee for bringing this fact to our attention. 



mf ) (x 1 ,...,x n )=Ef xij Xn) 
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following probability measure on the space of real N x N matrices: 



fi{dM) 



1 



lT,(M-M) Aiv><jv(dM) 



(3) 



(2tt) n2 / 2 



e 2 



where XnxN is Lebesgue measure on R . Even though the real Ginibre ensemble is 
a classical matrix model, the eigenvalue correlation functions have been computed only 
recently, see [5], [H], [9], [E] and [[6]. 

It turns out that the pfaffian point process corresponding to one-dimensional distri- 
butions of annihilating Brownian motions is equivalent to the pfafnan point process 
describing the law of real eigenvalues of Ginibre in the limit N — > oo. Namely, com- 
paring the statement of Theorem [2] with Corollary 9 of [5] we arrive at the following 
conclusion: 

Corollary 4 The one-dimensional law of particle positions for the system of annihilating 
Brownian motions on M at time t > under the maximal entrance law is a Pfaffian point 
process with the kernel 



where K^ r mtbre is the N — > oo limit of the kernel of the Pfaffian point process character- 
izing the law of real eigenvalues in the real Ginibre ensemble. 

In other words, the one-dimensional law of annihilating Brownian motions under the 
maximal entrance law initial conditions is equivalent to the N = oo limiting law of real 
eigenvalues of a real matrix with independent normal entries. 

Corollary[5]suggests that real eigenvalues of real matrix-valued Brownian motion might 
behave like a system of one-dimensional annihilating Brownian motions. In fact, numer- 
ical evidence we accumulated up to date compels us to make the following conjecture. 

Conjecture. Under the maximal entrance law, all finite- dimensional distributions of 
particle positions for a system of annihilating Brownian motions on R coincide with 
N — > oo limit of multi-time correlation functions of real eigenvalues of gl-R(N) -valued 
Brownian motion. 

Here gl-R(N) denotes the linear space of all N x N matrices with real entries. 

2 Brief review of some facts for one-dimensional coalescing 
and annihilating Brownian motions 

We consider, at first, initial conditions that have only finitely many particles. This paper 
describes only the one dimensional time distributions, that is at a fixed t > 0, of any 
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remaining particles. We list the positions of the particles at time t as (X\ : i E It). 
The exact details of the labeling system It will not be important for us, and indeed our 
results all relate only to the empirical measure N t defined by 

N t (A) = ^2 xi x t G A) f° r measurable iCR. 

For the case of annihilating particles, if the initial number of particles is even then it 
remains so for all time. To a list (xj) of an even number 2n of disjoint positions we may 
associate the open set 

S{(Xi)) = (xi,X 2 ) U . . . U (x 2 n~l,X 2 n) 

where x% < . . . < xm are the ordered positions. Some of the formulae for annihilating 
particles are then most easily stated in terms of the set valued process 

S t = S{(xi-.ieI t )). 

Notation. We write PP x ^ to indicate that we are considering (instantly) coa- 
lescing Brownian motions started from initial positions xi, . . . , x n . When the particles 
are annihilating we change the superscript from C to A. When the initial positions are 
random we change the subscript to 3, where S is the law of (Xq : i 6 Iq). 

2.1 The thinning relation 

The formulae about coalescing systems in the paper will always come with an analogue 
for annihilating systems. The close link between the two systems has often been observed. 
For this paper the formulae can usually be derived from the following thinning relation. 
For a list of positions (x\, . . . , x n ) we let 0(xi, . . . , x n ) be the random subset formed by 
thinning at rate 1/2, that is by removing each position independently with probability 
1/2. We may also thin a random set of positions, for example ®(X\ : i £ It), with the 
understanding that the randomness in the thinning is independent of the randomness 
in the set of positions. We write 0(E) for the law of the thinned set of positions that 
initially have law H. Then the thinning relation between coalescing and annihilating 
Brownian motions is the following equality in distribution: 

(XI :i€l t ) under P A {S) = 0(X| : * € I t ) under P§. (5) 

Such a thinning relation is discussed in Arratia p] for the scaled limit of reacting random 
walks, and is related to results in many later papers. There is a simple colouring proof 
(see ben Avraham and Brunet [3]) that bears repetition here. After the paths of a coa- 
lescing system have been realized, independently add random colours as follows. Initially 
colour each particle red or blue independently with probability 1/2. At coalescences the 
colours evolve according to the rules R + R — > R, B + B — )■ R and R + B — > B. Then 
the resulting system of blue particles evolves as an annihilating system. Moreover the 
colour of a particle at time t depends on whether there were a odd or even number of 
ancestors at time zero that were coloured blue. Since distinct particles have disjoint sets 
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of ancestors, the colour of all particles at any time t > remains independently red or 
blue with equal probability. The thinning relation follows. This argument makes it clear 
that the result holds much more widely, since the exact nature of the motion process is 
not relevant, nor is the mechanism of reaction (for example it holds for delayed reactions, 
when the reactions are controlled by the intersection local time). 



2.2 Duality formulae 

We use two duality formulae. For at < a 2 < ■ ■ ■ < a^m let Ik = (afc,afc+i) for k = 
1, . . . , 2m — 1. Then for disjoint (xj) 

P £) Wtih) = N t (h) = ■■■ = N t (I 2m ^) = o] = p ( i i} [S t n (xi) = 0] . (6) 
The annihilating analogue of this is, writing \A\ for the cardinality of a set A, 

~(_l) Ar t( / i u/ 3U...u/2m-i)" _ fiA "(_]_yStn(xi)f _ (7) 



There are various ways to see these formulae, but for coalescing systems a key construc- 
tion is the Brownian web and its coupling with the dual Brownian web, first considered 
by Arratia and explored in Toth and Werner |20j (and subsequent papers). We need only 
part of the Brownian web as follows. For a fixed t > 0, there is a system of coalescing 
Brownian motions starting from every rational x and running over the time interval [0, t], 
and a coupled system of backwards coalescing Brownian paths starting at time t at all 
x G Q and running back to time zero. In fact, the Brownian web has particles starting 
at all space-time points (s, x) but we will not need this, and it is enough to establish 
([6]) first for rational (xj) and (cij). The key property is that, almost surely, none of the 
forward paths cross any of the backwards paths. (A discrete version of this coupling, that 
is using simple coalescing simple random walks, is easy to construct - see the appendix 
in [18] - and illustrates this non-crossing property). From this non-crossing property one 
sees that the event that N t ((a, b)) = for the forward coalescing system is almost surely 
equal to the event that the open interval formed by pair of backwards particles starting 
at a and b does not contain any of the initial forwards particles. The coalescing duality 
([6]) follows immediately, once one notes that St may be replaced by its closure and that 
annihilating the backwards particles when they meet will not affect this closure. 



The annihilating duality ([7]) follows from ([6]) and the thinning relation. Note that 
thinning a set of n > 1 elements produces a random subset whose size has a binomial 
B(n, 1/2) distribution, and also that E[(-l) B ^ n ' 1 / 2 ^} = 0. Then thinning and © show 
that 



E, 



e(xi) 



_ 1 \A r t (/lU/ 3 U...U/ 2 m-l) 



P(t ) [s t n(x i ) = ®] 



N t (L 



2m- 1 j 



0] 



E 



(Oi) 



.-^IStneooi 



(where on the right hand side ^ is the expectation over the annihilating particle 
system and over the independent thinning). One may then argue by induction on the 
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number n of the initial particles (xi, . . . ,x n ). When n = 1 the above identity reduces 
to ([7]) for a single particle. For general n the identity is a mixture of copies of ([7]) for 
initial conditions that are subsets of (x,). But all but one of the copies will involve n — 1 
or less particles allowing an inductive proof. Note also that ([6]) also follows from ([7]) - a 
weighted sum of dTJ according to the distribution of 0(a?j) yields (0). 

Remark. Other coalescing duality formulae, such as those in Xiong and Zhou |21j . 
also follow from the Brownian web and its dual, but their proof shows that one may also 
establish them using the Markov generator duality, as explained in section 4.4 of Ethier 
and Kurtz, and thus bypass the Brownian web. In particular this generator technique 
may be extended to show analogous dualities for more general spatial motions, where 
the web construction does not (as yet) exist. Formally the generator proof shows that 
the dualities ([6]) and © will hold for instantly reacting continuous Markovian motions, 
where the motion on the right hand side must be the image of the motion on the left 
hand side under reflection x — > —x. Furthermore the maximal entrance laws constructed 
in the next section should follow once some moment control is established, which will 
require some non-degeneracy of the spatial motion to ensure enough reactions take place. 

2.3 Maximal entrance laws 

One may start coalescing systems from infinitely many particles at time zero. A natural 
state space for the empirical measure is the set M lfp(R) of locally finite point measures 
on R, which is a closed subset of the space of locally finite measures under the topology 
of vague convergence of measures. The reactions ensure that the point masses only have 
mass one, and so we consider the (measurable) subset A^o of those measures of the form 

fj, = ^^S Xi where (xi) is locally finite in R and has disjoint elements. 

i 

(To obtain a process with continuous paths, which does not concern us in this paper, 
one can quotient Mlfp by the minimal relation that ensures fi + 25 x ~ fj, + 5 X .) 

There is a Feller Markov transition kernel ptif^jdv) on Mq. Moreover, there is a 
maximal entrance law, intuitively starting with one particle at every site (as in the 
Brownian web). This can be characterized by passing to the limit in ([6]) as (xi) increase 
to become dense in the real line. This entrance law, which we denote by Pg, has one 
dimensional distributions satisfying 

Pg [N t (h) = N t (I 3 ) = ... = iVi(Wi) = 0] = Pfa [t < t] (8) 

where r is the time for complete extinction of the annihilating system. This characterizes 
the one dimensional laws on .Mo, and these laws are an entrance law for the Markov 
transition kernel described above. By the scaling property of Brownian motions we have 
P(Ta-)l T ^ T 2 t] is independent of T > 0. Using ([8]) this translates into a scaling for the 
entrance law 

The law of (T^Xj 2 : i £ I tT2 ) is independent of T > under Pg. (9) 
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Many suitably spread out and non-degenerate initial conditions are attracted to the 
maximal entrance law as t — > oo. For a large class of initial conditions (xj), the law of 
(T~ 1 X l T2t : i € I^t) under PEa converges in distribution, on Mlfp(R) as T — )• oo, to 
the law of (Xf : i G It) under P~,. Indeed, using the extension of ([6]) to countable (xj), 
this follows (see the appendix) from 

Pg.) [(T- 1 X^ t : i 6 I T 2 t ) fl J fc = for k = 1, 3, . . . , 2m - l] 

= Pg [A^CT/i) = 7V tT2 (T/ 3 ) = • . • = N tT 2(Tl2 m -x) = 0] 

= p ( 4, aj) [5 tT2 n (xi) = 0] 
= p ( t) n (rt.) = 0] 

= P^ [JV t (Ji) = iV t (J 3 ) = . . . = N t {hm-i) = 0] . (10) 

The third equality comes from Brownian scaling and the final equality is ([5]). The 
convergence holds for deterministic (xj) for which (T _1 Xj) become dense in any finite 
interval [a, b] as T — > oo. A large class of random initial conditions will clearly also work, 
for example non-zero stationary and spatially ergodic. 

For annihilating systems a Markov transition kernel can also be constructed, using ([7]) 
and it's extension to countable (xj) as a means of characterization. We can define an 
entrance law P^ for the annihilating system by taking the thinned copy of the entrance 
law for the coalescing system. This satisfies the formula 

= P(i) [r < t] (11) 

which again determines one dimensional laws on Mq that form an entrance law for the 
annihilating system. The domain of attraction of this entrance law is more delicate. The 
example in section 3 of Bramson and Griffeath [6] suggests that different approximations 
to a maximal entrance law may yield different laws at times t > (their example uses 
varying intensities of nearby pairs at time zero). For initial conditions that fill the 
lattice A _1 Z, or that are Poisson with intensity A, the one-dimensional time distributions 
converge as A — > oo to those of the entrance measure, or for a fixed A the large time 
distribution rescales to those of the entrance law, by the argument above. 

Since we found it difficult to find a full account in the literature, we give, in the 
appendix, a brief sketch of the proofs of the results in this subsection. 



3 Proofs 

3.1 Review of PfafRans 

We give a short summary, targeted at beginners like us, of the facts we shall use about 
Pfaffians (mostly proved in [19] section 2), and of the definition of a Pfaffian point 
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process. We write Pf(a y - : 1 < i < j < 2n) (or just Pf (a^ : i < j)) for the Pfaffian of 
the real anti-symmetric matrix whose elements are for i < j. 

The determinant of an anti-symmetric matrix of odd order is zero. Suppose A is an 
anti-symmetric 2n x 2n matrix. Then det(A) is the square of a polynomial of degree n in 
the matrix elements, called the Pfafnan of A and written as Pf (A). One can define the 
Pfafnan as a suitable sum over permutations of products of matrix elements. Indeed, 



(12) 



where T<2n is the set of permutations a of {1, 2, ... , 2n} given by a{2k — 1) = ik, <r{2k) = 
jk for k = 1, . . . , n for which the choices (ifc), (j^) satisfy i& < for all k and i\ < 12 < 
. . . < i n . A convenient way to calculate the sign of such a permutation is via crossings. 
The quadruple ik,jk,ii,ji is called crossed if ik < %i < jk < ji- Then the sign of a 6 E2 n 
equals (— 1) M where M is the number of crossings. To visualize these crossings easily 
one can embed the integers 1, . . . ,2n into the x-axis of the plane and join ij~ to jk for 
each k with a loop in the upper half plane. 



It is worth recording the smallest cases: 

Pf \ ' '.. =af-be + cd. 

The explicit 4x4 case was used to guess many of the Pfafnan formulae in this paper. 



Pf 



a 
-a 



= a 



( a b c \ 

—a Ode 

-b -d / 

V -c -e -/ J 



Pfaffians have many similar properties to determinants. It follows from the definition 
that Pf(AjAjajj) = Pf(a y ) n^-^fc- Por an y ^ n x ^ n ma tri x B the product B T AB is 
anti-symmetric and Pf (B T AB) = det(S)Pf (A). The Pfaffian can be decomposed along 
a row, or column, of the matrix. For example if A is a 2n x 2n anti-symmetric matrix it 
satisfies the recursion, for any i € {1, 2, ... , 2n}, 



2n 



Pf (A) = {-l) i+j+1 aijVi{A^) 



(13) 



where A^'^ is the (2n — 2) x (2n — 2) submatrix formed by removing the ith and jth 
rows and columns. We will also use a decomposition formula for the Pfaffian of a sum 
of two 2n x 2n anti-symmetric matrices A and B, namely 



Pf (A + B) = J](-l)l J l/ 2 (-l) s ( J )Pf (A\ j)Pf (B\ jc) 



(14) 



where: the sum is over all subsets J C {1,2, . . . , 2n} with an even number of terms; 
J c = {1,2, . . . ,2n} \ J; s(J) = ^2j e jj (and s(0) = 0); and where A\j means the 
submatrix of A formed by the rows and columns indexed by elements of J (and the 
Pfaffian of the empty matrix is taken to have value 1). 
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Suppose a measurable kernel 



K[X ' y) V K 21 {x,y) K 22 (x,y) J 

is anti-symmetric, in the sense Kij(x,y) = —Kji(y,x) for all i,j £ {1,2} and x,y £ 
R. Suppose it also acts as a kernel for a a bounded operator on L 2 (R) © L 2 (R). A 
point process (X 1 : i & I) with n-point density functions p <yn \x\, . . . ,x n ) is called (see 
Soshnikov |17j ) a Pfaffian point process with kernel K if p( n \x\, . . . , x n ) is given by the 
Pfaffian of the 2n x 2n anti-symmetric matrix formed by the n 2 two-by-two matrix entries 
(K(xi,xj) : i, j = 1, . . . , n). The kernel is not uniquely determined. 

A very convenient tool for manipulating Pfaffians is the Berezin integral. We provide 
arguments that avoid this tool in this paper, and so do not describe the rules for manip- 
ulating these integrals. However they were used repeatedly while exploring these results, 
and in the next section we show how the Berezin integral can considerably shorten the 
argument. A very readable account of Berezin integrals can be found in Itzykson and 
Drouffe |llj . The key property linking the Berezin integral to Pfaffians is (compare with 
the normalizing determinant for multi-dimensional Gaussian integrals) 

Pf (ciij :i<j) = J d^ 2n ■ ■ ■ #ie"5 £&=i ^ a «^ . (15) 



3.2 Proof of Theorem 02, the product moment Pfaffians 



We start with the product moment, defined for bounded measurable g : R — > R and 
disjoint (xj) by 

ieh 

where the product over an empty set, occurring when all the particles have been annihi- 
lated, is defined to have value 1. Note that m^(x) is given by the Brownian semigroup 
applied to g. We also set = 1. To show that mS 2 ^ is given by a Pfaffian, we shall 
give a p.d.e. derivation similar in spirit to that showing the Karlin McGregor formula for 
the transition density for non-intersecting Brownian motions is given by a determinant. 

Let V n C R n be the open cell {x : x\ < x 2 < ■ ■ ■ < x n }. On (0, oo) x V n the function 
m[ n \x) solves the heat equation, and we must examine the boundary conditions. For 
n > 2 and when g is bounded and continuous, the functions are continuous on 

[0, oo) x V n and extend to a continuous function in C((0,oo) x V n ). There are lots 
of pieces to the boundary of V n , but the most important are the faces Fi^ n defined by 
Xi = Xj+i and where the remaining x^ are disjoint. On F{ n the continuous extension 
agrees with the lower order moment m*- n_2 -'(a;^' i+1 )), where x™' £ R n_2 is the (n — 2)- 
tuple formed by removing X{ and Xj from (x%, . . . , x n ). This can be seen by showing that 
near the boundary the hitting time between particles starting at Xi and Xj+i is likely to 
occur before any other collision and before time t. On other parts of the boundary the 
extension agrees with other lower moments. 



m\ n \x 1 ,...,x n ) =Ef Xi> ^ Xn) 
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The system of heat equations for (m^ : n = 1, 2, . . .) 

^m ( j n \x) = Amf'(x) on (0, oo) x V n , 

m[ n \x) = m[ n 2 ^(x^ i,i+1 ^) for x € and i = 1, . . . , n — 1, 
m o n) ( x ) = nr=i5'(^) for x 6 V ni 

forms a closed system, in that each equation has boundary conditions formed by equa- 
tions of lower order. Note that, typically, the initial condition does not match the 
boundary conditions. Taking g bounded and smooth, the system has unique solutions 
in C 1,2 ([0,oo) x V n ) n C((0, oo) x V n ). It is enough to specify boundary conditions only 
on each face F^ n - the Feynman-Kac formula makes it clear that the other parts of the 
boundary of V n do not affect the value of mS n >. 

To establish the Pfaffian ([2]) stated in Theorem [3j it is enough, by an approximation 
argument, to treat the case where g is smooth. We shall prove ([2]) by showing the Pfaffian 
Pf (m^ 2 \xi,Xj) : 1 < % < j < 2n^ solves the system of heat equations above. Note that 
(|2j) holds when t = since 

2n 

mf n \xi,. . . ,x 2n ) = ~[]_g(xi) = Pf (g{xi)g(xj) : i < j) . 
i=i 

The Pfaffian is a finite sum of product terms (see (112p ) of the form 

sgn(o-)mf } (x h , x h )mf ] {x l2 ,x h )... mf ] (x ln , x jn ) 
where a is a permutation given by a (2k — 1) = if., a(2k) = jj. for k = 1, . . . , n. Since 

(2) 

m\ (x,y) satisfies the heat equation on [0, oo) x {x < y}, each product term lies in 
C 1,2 ([0, oo) x V^n) and satisfies the heat equation on [0,oo) x y 2n . Since m t (x,y) 
extends continuously to (0, oo) x {(x,y) : x < y}, the Pfaffian extends continuously to 
(0, oo) x V2n- By the uniqueness for the system of heat equations, it remains to check 
that the Pfaffian satisfies the required boundary conditions on each face Fi^n which will 
complete the proof of Theorem [3j 

We show the argument for the face F\^ n where x\ = x 2 (other faces are similar). 
We may argue inductively, and suppose that m^ k ' is given by the Pfaffian for k = 
0, 2, . . . , 2n — 2. Our quickest proof is using the representation (|15|) in terms of Berezin 
integrals. This gives 



Xl=X 2 



X\=X 2 



Pf (mf\xi,Xj) : i < j 

. . . ^ ie -|E^3^ m P ) (^^)^ e -Wi+^)E'= 3 ^ > (^i^fc)^. 



The sum M = Y^kLs r ni 2 \ x i' x k)' l Pk is independent of ipi and ip 2 and the dip 2 dip\ integral 
becomes (using the rules for Berezin integrals) 

j d^ 2 d^e~^ + ^ M = J # 2 #i(l - V^iXl - (V'l + fh)M) = 1. 
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This simplification leaves only J dip2n ■ ■ ■ dtp^e 2 S»,j=3^i m t (?H&j)il>j which [ s the Berezin 
integral for m( 2n_2 )(x3, . . . , X2n)- 

An argument that avoids Berezin integrals is as follows. Using the recursive relation 
for Pfaffians (|13p we see that the Pfaffian in (|2|) equals 

2)i 



y~^(— l) k mf\xi,Xk)mf !./- 



k=2 



Since m\ (x\, X2) extends to the function 1 on x\ = X2, it remains only to check that 

2n 

^(-l^mf^xx^mf^V 1 '^) (16) 

fc=3 

vanishes when xi = X2 and i > 0. But this follows from expressing m^ 2n ~ 2 ^ using (|12p . 
Indeed, fix j, k > 3. Then for an expression of the form 

(2)/ \ (2)/ \ (2)/ s (2), \ 

; (xi,x fc )m^ '(afc.ajj-Jmi (x i2 , x j2 ) . . . ; (x in _ 1 ,x jn _ 1 ) 

arising from the feth term in (fl~6j) . where {^2,^2, • • • > *n-i, Jn-i} = {3, 4, . . . , 2n} \ {j, fc}, 
there is a corresponding term 

(2)/ \ (2)/ s (2)/ v (2)/ \ 

ml '{xx^jjml >{x2,x k )m\ '{x l2 ,x h ) . . .m\ '{x^, x jn _J 

arising from the jth term in (|16p . These terms agree on x\ = X2 and a careful check of 
the signs of the permutations, and the factors (— 1) J and (— l) k in (|16fl . shows they will 
cancel. One way to do this check is to compare the sign of 

_ / 2 3 4 5 ... k - 1 fc + 1 ... 2n - 1 2n 

\2 j ! 2 !3 in-1 in-1 

with that of 

, ( 2 3 4 5 ... j-1 j + 1 ... 2n - 1 2n 



2 k %2 h in-l jn-1 

by counting crossings. The loop joining 2 to j in a must be replaced by a loop joining 2 
to k in <r'. This may affect crossings with any of the loops emanating from sites between 
j and k, and will do so unless a pair of them are joined to each other. There are \k— j\ — 1 
sites between j and k so it will change the parity of the number of crossings exactly when 
\k — j\ is even. 

Remark 1. For odd moments there is also a Pfaffian representation, namely, when 

Xi < X 2 < • • • < X2n-1, 

(2n-l), m \_v>rf^ (2) 



(xi, . . . ,x 2n -i) = Pf (mfix^Xj) : < i < j < 2n- lj (17) 

where we adopt the convention that (xQ,Xk) = m\ (x^). This Pfaffian involves a 
linear combination of terms of the form 

sgn(cr)m^ 1) (x jl )m (2) (x i2 ,x i2 ) . . . m {2) (x in , x jn ) 
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which again shows that it solves the heat equation when [0, oo) x Vm—\- The recursive 
Pfaffian relation gives 

2n-l 
k=l 

Expanding the Pfaffian along its first row using (fT3"j) we obtain for x = (xi, . . . , X2 n -i) £ 

2n-l 

Pf (mf\x uXj ) : < % < j < In- l) = {-l) k+1 mf ] {x k )mf n ~ 2 \x^) (18) 

fc=i 

where we again write superscripts x^'^'"' to mean that we remove the indicated co- 
ordinates. The terms with k = 1 and k = 2 cancel on the face i*2n— 1 1 where x\ = X2- 
Moreover on this face, for k > 3, m[ 2n 2 \x^) = m[ 2n ^ (a^ 1 ' 2 '^)) so that the Pfaffian 
in (fT51) becomes 

2n-l 

J2(-l) k+1 rn^{x k )ml 2n - A \x^). 

k=3 

But this is the decomposition of va\ n 3 \ X ( 1,2 ^) along the first row, and this shows the 
boundary conditions are correct on i*2n— Other faces are similar. 

Remark 2. Since our proof relies only on uniqueness for the underlying system of 
heat equations, the extension of these product moment Pfaffians to more general spatial 
motions looks quite straightforward, for example to Markovian spatial motions that are 
continuous and suitably non-degenerate. The Pfaffians in the next section would then 
also follow for these more general motions, just by algebraic manipulation, once maximal 
entrance laws characterized by © and (jlip are established. 



3.3 Proof of Theorem [2, the Pfaffian point process kernel 



Fixing di < ... < a2m and choosing g(x) = (— l)S<x(a ; <ai) [ n (JJJ) we see that both sides 
of the duality ([7]) are Pfaffians in the variables (xi). Choosing g = 0, recalling that an 
empty product takes the value 1, we see that Pu..\[ T < is a Pfaffian. The entrance law 



dualities ([8]) and (fll) show that 

pg[N t (h) = N t {h) = . . . = N t {L 



2m- 1, 



0}=E t 



\N t (hUI 3 U...UI 2m - 1 ) 



are Pfaffians in the variables (a^). The entries in this last Pfaffian are explicit since 



Pg[N t ((a J ,a k )) = 0] 



Ei 



1) 



N t ((aj,a k )) 



are all equal to (by Brownian hitting probabilities) 

F (t _1/2 (aj - a;)) where F(x) = tt' 1 / 2 J™ exp(-y 2 /4) dy. 



(19) 
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We switch dummy variables for the rest of this section, taking x\ < X2 < ■ ■ ■ < £2n and 
Ik = (%k,Xk+i) so that we start from 

Pg [N t (I k ) = for k = 1, 3, . . . , 2n - 1] = Pf (f (^ 1/2 (^- - x f )) : i < j) . (20) 

To prove Theorem [2j we shall identify the Pfaffian point process kernel by differentiat- 
ing the empty interval Pfaffian (|2U|) above. By scaling we may take t = 1. Differentiate 
the identity (|20p for t = 1 in the variables x±, X3, . . . , x 2 n— 1> The left hand side becomes, 
formally, 

[iVx(dxi) JVi(dx 3 ) • • • iVi(dx 2ri -i)I(^i(4) = for k = 1,3, . . . , 2n - 1)] . 

Letting x 2 i | X21-1 for Z = 1, . . . , n we reach the n-point density pf'\x\, x$, . . . , a^n-i)- 
In the appendix 14,31 we give more details verifying the formal differentiation above is 
valid, by using distributional derivatives. 

On the right hand side of (|20p we will also differentiate in the variables 2:1,2:3, ... , x 2n -i- 
Note that each product term in the permutation expansion (112j) of the Pfaffian contains 
exactly one element that involves the variable x\. So differentiating in X\ leads to a 
similar permutation expansion, but where all the terms that involve x\ have been differ- 
entiated. Repeating this argument, differentiating in x\, a; 3 , . . . , x 2n -\ yields the Pfaffian 
where each term in the matrix has been differentiated in the variables £1,2:3, ... , X2n-i 
that is where the 2x2 block formed by the rows 2j — 1, 2j and columns 2k — 1, 2k is 
given by 

/ -F"(x 2 k-1 - %2j-l) -F'(x 2 k - %2j-l) \ 

V F'(x2k-i ~ X 2 j) sgn(x 2 fc - X 2 j) F{x 2 k - X 2 j) J 

when j < k. (Note that F" is an odd function and so no sgn is needed in the 2j — 1, 2k— 1 
entry.) Letting X21 i X21-1 for I = 1, . . . , n we obtain the kernel K stated in Theorem [2j 
The decay in F, F' , F" implies that K acts as a suitable bounded operator. The scaling 
relation (|9|) implies that the kernel of the distribution time t is t~ 1 ^ 2 K(xt~ 1 ^ 2 , yt^ 1 ^ 2 ). 

Remark 1. An alternative starting point, used by ben Avraham et al. (see [2], [3]), 
is to show the empty interval probabilities Pg[N t (Ii) = Nt(Is) = . . . = N t (l2m-i) = 0] 
satisfy heat equations in the variables (x{), though the connection with Pfaffians does 
not seem to have been noted. 

Remark 2. The linear ordering of particles seems to be crucial. For Brownian particles 
on a one-dimensional torus, there is an extra boundary condition where X2n may hit x\ 
by going 'the other way' around the torus, and this is not satisfied by the Pfaffians. 

Remark 3. Differentiating in the variables x 2 , ... ,X2 n instead leads to the alternative 
dual kernel 

/ sgn(y - x)F(\y - x\) F'(y-x) \ 
{ ,V) \ -F'(y-x) -F"(y-x) J 

Remark 4. Starting from the Pfaffian (|20p certain other probabilities can, by algebraic 
manipulation, also be expressed as Pfaffians. We give three examples, leaving details of 
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the derivations to the appendix. In each case F is the 2n x 2n anti-symmetric matrix 
with elements Fij = Pg[N t ((xj,x k )) = 0] = F(t~ 1 / 2 (xj - Xi)) as in ([15]) . 



Let I = l2n be the 2n x 2n anti-symmetric matrix with entries 1 above the diagonal. 
Then 

Pg [N t (I k ) > for k = 1, . . . ,2n- 1] = Pf (/ - F). (21) 
The annihilating analogue of this is 

P* [N t {I k ) is odd for k = 1,2, . . . ,2n - 1] = 2 1 ~ 2n Pf (/ - F). (22) 



matrix ( ^ ) down the diagonal and zeros elsewhere. Then 



Let O = 02n be the 2n x 2n anti-symmetric matrix formed by n copies of the 2x2 
1 

1 

Pg [N t (h)) > for k = 1,3,5,... ,2n- 1] = Pf (O - F). (23) 
Again there is an annihilating analogue. 
• Let O = 0<m be the 2n x 2n anti-symmetric matrix with entries 

{+1 if i = 2,4, . . . ,2n - 2 and j = i + 1, 
-1 if j = 2,4,...,2n-2andi = i + l, 
otherwise. 

Note that O also has copies of the 2x2 matrix y ^ ^ ^ in some places down 
the diagonal and zeros elsewhere. Then 

Pg [N t (I k ) = for k = 1, 3, . . . , 2k - 1 and iV t (4) > for k = 2, 4, . . . , 2k - 2] 
= Pf(F-O). (24) 

Remark 5. Suppose (M xy : x < y) is a bounded continuous field that satisfies 
^[M^^M^J = Pf (£/[Ma. iXj ] : 1 < i < j < 4) for x x < . . . < x A . Then by continuity 
E[M XX M XX ] = Pf (E[M XX ]) = E[M XX ] 2 and so M xx must be deterministic. This imposes 
a restriction on the class of correlation functions admitting a PfafRan representation. 

3.4 Proof of Theorem [TJ the asymptotics for pf n ^ 

We work throughout under the entrance measure Pg. By thinning the corresponding 
density for annihilating systems differs only by a multiplicative factor 2~ n . The n-point 
density function p[ n \x) is a Lebesgue density for the measure E'glN^dxi) . . . N t (dx n )} 
on V n . The existence of this density, defined almost everywhere, and the simple bound 

Pt n \x) < C n t' n/2 for alH > and x € V n (25) 
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is discussed in [12]. Furthermore there we established the following upper bound: for all 
L > there exists Cl < oo so that 

p\ j (xi,...,x n ) < C L t~2 3— J | \xi-Xj\ for all t > and \xi\ < Lt 1 ' 2 . (26) 

l<i<j<n 

As t — > oo the entries in the Pfaffian for are of the form F, F', F" evaluated at 
points t~ l l 2 (xj — x{) close to zero. One may approximate these by using the Taylor 
expansion for F(z) at small values of z. However, considerable cancellation occurs in the 
many terms of the Pfaffian and it is not immediately clear how to read off the leading 
asymptotic decay in t. Indeed the following argument shows at F needs to be expanded 
to a large number of terms to obtain the correct answer. 

We shall analyze first a modified density function pf n \x) for x £ V2 n , which is a 
density for the measure 

Eg [N t (d Xl ) . . . N t {dx 2n ) x (N t (I k ) = for k = 1, 3, . . . , 2n - 1)] 

(where we recall that 1^ — (%ki ^fc+i 

)). We claim that 

~pf n \x u x 2n ) = (47rt 2 )-"/ 2 Pf (V ({ Xj - x 4 )/t 1/2 ) : 1 < i < j < 2n) (27) 

where 4>{z) = z exp(— z 2 /A). This follows formally, as in section T3.31 by differentiating 
(f20|) in all variables x\,x 2 , ■ ■ ■ , X2 n , and using that, 

d Xl ... d X2n Pf {F( Xj - Xi ) :i< j) = Pf (-(4vr)- 1 / 2 0(x i - a*)) 

(which follows from differentiating each term in the permutation expansion (| 12j) of the 
Pfaffian). We give more details in the appendix 14.31 



The advantage of the representation (|27p is that it is a Pfaffian all of whose entries are 
of the form /(xj — Xj) for a single function /, and this allows us to apply the following 
lemma, proved at the end of this section, that gives an expansion for a Pfaffian whose 
entries are close to the zero of an odd function. 

Lemma 5 Let (j) : R — > R be an odd function that is analytic at zero. Then for any 
n > 1 there exist e(n, <fi) > and C(n, 4>) < oo so that for y £ V n with \y\ < e(n, cp) 

Pf (HVj -Vi):l<i<j<2n) = Pf(J {2n) (<P) + R {2n) (y)) U (Vj ~ 2/0 

l<*<j<2n 

where J^ 2n \(p) is the constant anti- symmetric matrix with entries 

rW'Hr' ^^.^^ IO) forl<i< j <2n, (28) 

and the remainder R^ 2n ^ (y) is a anti- symmetric matrix satisfying 

R^ n \y) < C(n,4>)\y\ for all i,j and \y\ < e(n,</>). 
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We apply this lemma to the Pfafflan in (|27p with <p(z) = ze z2//4 and with for 
t large enough. Expanding the Pfaffian Pf (j( 2n )(0) + J R( 2n )(t~ 1 / 2 x)) using (JT3|) we find 
only one term, namely Pf (j( 2n )(0)), that does not decay as t — > oo. This shows that 

lim e 2 ^pf n \x u ...,x 2n ) = (4vr)~"/ 2 Pf (j( 2? %)) TT ( Xj - Xi ). 

l<i<j<2n 

To obtain the same estimate for p( 2n ) we estimate the difference as follows. 

< p[ 2n \xi, . . . ,X 2n ) ~ pf n \xi, . . . ,X 2n ) 

= Eg \N t {dxi) . . . N t {dx2n)x(N t {Ik) > for some k = 1, 3, . . . , 2n - 1)] 
2n-l 

< E ™ [Ntidxi) . . . N t (dx 2n )N t (I k )} 

k=l 
2n-l 



/ Pt (^i,---,Sfe,^,^fc+i,---,a;2n)^- 

fc=i 



Each term in this sum is of a smaller order in t by 

Examination of the proof shows that we need not let the values of x\, . . . , X2 n be fixed, 
and that in fact we may take the supremum over any positions (xi(t)) provided that 
supj |xj(i)|i -1 / 2 — > as t — > oo. 

Proof of Lemma [5], Let $ be the 2n x 2n anti-symmetric matrix with entries given 
by &ij = 4>{yj — yi). The aim is to show, for small y, that 

$ = V T (J + R)V 

where J and R are as in the lemma (with n fixed and suppressed) and V is the 2n x 2n 
Vandermond matrix given by V!y = Vj~ ■ Since det(V) = Yii<i<j<2n(yj ~ Vi)> the 
conclusion then holds from Pf (V T {J + R)V) = det(V)Pf (J + R) " 

For small \y\ we expand by analyticity (writing (j) k (0) for the fcth derivative of 4> at 
zero) 



oo 



n=0 



oo n 1 



n=0 fc=0 



oo 

E u7f^(o)^(-^)' 



fc,Z=0 
oo 

fc,Z=l 

where we have rearranged using / = n — in the penultimate equality. Note that 

2n 2n 

fej=i fej=i 
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It remains to re-express the remaining terms in (|29p as the desired remainder. 
Recall the symmetric polynomials cr^ n (y) defined for y £ R 2n by 

2n 2n 

H(y k - A) = ^(-l)Vr(y)A 2 «- fc . (30) 

k=l k=0 

Note that a^ 1 is a polynomial of order k. Since a^ 1 = 1 we may choose A = yi to see 
that 

2n 

= vT + ^(-l) fe af (y)j/f - fe for i = 1, . . . , 2n. 

k=l 

Multiplying by yf we see that 

2n 



y 



Y^(-l) k+1 a 2 k n (y)yf +2n - k for i = 1, . . . , 2n and p = 0, 1, . . .. (31) 



fc=i 



By iterating this we may express yf +2n for p > as a mixture of 1 , y% , yf , . . . , y 2n 1 , as 
follows: 



2n 

2/f +2n = £ ^^'"(y)^" 1 far i = 1, . . . , 2n and p = 0, 1, . . . (32) 
fc=i 

where T^ n ' p+2n (y) is a polynomial of order p + 2n — k + 1. Using this substitution in the 
remaining terms of (|29|) . that is where k or / is at least 2n + 1, we find (formally) that 

oo 2n oo 2n oo \ 2n 

£ +E E +E E = E f 1 ^) 

k fc,Z=2n+l fc=lZ=2n+l 1=1 k=2n+l J p,9=l 

where 

oo 

«w(v) = E r^ l -\y)J kl r^ k -\y) 

k,l=2n+l 

2n oo 2n oo 

+ E E ^- 1 (y)^ + EE^" 1 (y)V (33) 

fc=H=2n+l 1=1 k=l 

Note the lowest order of the polynomial entries in the terms for R pq is of order 1 . In the 
appendix 14.41 we check that this rearrangement of (|29p is valid when \y\ is suitably small 
and that the required error bound \R pq (y)\ < C(n,c/>)\y\ holds. 

Thanks. We would like to thank our colleague Dmitriy Rumynin for advice on the use 
of symmetric polynomials. 
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4 Appendix 



4.1 Details for section [2731 

We give a few details on (one approach to) the results surveyed in section 12.31 For 
coalescing systems one can use monotonicity, adding initial particles one by one, to 
construct the infinite system. This is not available for annihilating systems, so we sketch 
a weak convergence argument that applies to both systems. 

One can control moments by bounds on the n-point density function. Indeed p[ n \x), 
the density for the measure E^ x ^[N t (dxi) . . . N t (dx n )} on V n , depends on the initial 

condition, but satisfies the bound p[ n \x) < C n t~ n / 2 uniformly over all possible finite 
initial conditions (xi). This follows by duality for n = 1 and by anti-correlation for n > 1 
(see [32])- It follows that E? x ^ [Nf(a, b)] is bounded, for each t,p > 0, a, b G R, uniformly 
over finite initial conditions (xj). 

Fix \x G M-o and take finite measures (j, n so that p, n — > p (recall we are using vague 
convergence). The moment bounds above imply that the laws of Nt on Mlfp under 
P^ n are tight. Take a subsequence n' along which they converge to a limit, which we 
denote Q. The functions 

v -> F M (u) := xWi) = Hh) = ■■■ = Hhn-i) = 0) 

are discontinuous on M.lfp- However the moment bound E9 x n [Nf(a, b)] < C(t)(b — a) 
holds also for the limit law Q and implies that v({a,i}) = 0, Q{du) almost surely. This 
shows that Q does not charge the discontinuity set of Ft ai y Then we may pass to the 
limit in © to deduce that 

j F (ai) (v)Q(dv) = [S t n supp(^) = 0] . (34) 

These functionals do not characterize a law on M-lfp-, but they do characterize a law 
that is supported on Ado- To see this note that for u £ Mo 

*, ( k k + 1 \ 

v([x,y]) = jv lim w z^ x ( l/ ([ x ' y ] n (] v : '^^]) >0 ) • 

k ' 

From this one may use (f34"|) to find f v([xi,yi\) . . . v([x n , y n ])Q(dv) which, by the moment 
bounds, determine Q. To see that Q is supported on M.q note that 

P£ n [N t (a,b) > 2] < f ! pf\x 1 ,x 2 )dx 1 dx 2 < C(t)(b-af. 

J a J a 

This bound holds uniformly over n and hence also for the limit law Q. Then the conclu- 
sion follows from the usual covering argument, for instance 

Q [ p({x} > 1 for some x £ [— L, L\)\ 

LN 

< Q[^[k/N,(k + l)/N})>2\<C{L,t)N- 1 . 

k=-LN 
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Thus the law Q is determined and we may define pt(^,dv) to equal Q(dv). 

The remainder of the results in section I2TB1 follow using similar tools. For example, for 
the continuity of fj, — > pt(^, dv), that is the Feller property, suppose that [i n — > [x in .Mo- 
The moment bounds, which still hold for infinite initial conditions, imply the tightness 
of pt (/i n j dv) . Passing to the limit in 

J F{ ai ){v)Pt{Hn,dv) = Pfa [S t n supp(/i n ) = 0] . 

shows that any limit point of pt(p- n , dv) must be pt(^, dv). The semigroup property, for 
bounded continuous F : Mlfp — > R, 

j F(v)p t+S (n,dv) = J J F{v')p s {v,dv')p t {n,dv), (35) 

which is valid for finite measures /U extends to hold for [i E A4q by approximation, 
using the Feller property. The same tightness and characterization methods establish 
the existence of a law characterized by (jH|, and justify the arguments in (I10p that many 
initial laws are attracted to it. That (|8|) determines an entrance law can be established 
by passing to the limit in (j35|) along \i = ^ k <5 A -i fc as A — > oo. 

The annihilating case follows the same lines, with moments controlled since the n- 
point density and moments for annihilating systems are bounded by the corresponding 
coalescing system. The coalescing duality formula ([6]) is replaced by the annihilating 
duality formula ([7]), and to see that this will characterize the law note that for v € Mo 

v([x,y])= lim y (l-(-l)K[^]n(£,^])V 

N— >oo V / 
k 

4.2 Details for Remark 3 in section 13.31 

We give here the algebraic manipulations to derive the Pfaffians (I21l22l23l24p . 

Pf (/) = 1 for / the 2n x 2n anti-symmetric matrix with entries 1 above the diagonal, 
and the formula (fH|) specializes to 

Pf(I-A) = ^(-l) s ( J )pf(A| J ), 

(using for a 2n x 2n anti-symmetric matrix A, that Pf (—A) = (— l) n Pf (^4)). We combine 
this with a simple combinatorial identity (which can be checked by induction on n): 
suppose that (m^fc : 1 < j < k < n) satisfy the collapsing product rrij^mkj = rrijj for 
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all j,k,l; then 



n-l 



Y\ (1 + mjfc,fc+i) = 1+ ^2 m kiM 

k=l l<fci<fc2<n 

+ ^ rn klik2 m k3 ,k 4 

I<ki<k2<k3<k4<n 

+ ^ rn kl!k2 m k3:k4 m k5jk( . + ... 

I<ki<k2<k3<k4<k^<k^<n 

j 

where the final sum is over all subsets J of {1,2, ,,.,n} of even size, and if J = 
{ki, . . . , k 2m } where ki < . . . < k 2m then mj = m klik2 m k3:k4 . . . m k2m _ 1)k2m (and with 
m0 = 1). If n is even then the last term of this series is mi 2^3,4 . . . m n _x iTl . Note that 
m j,fc = & k ~ : * m j,k also satisfy fhj yk fh ky i = fhjj and applying the above for m one obtains 
a decomposition for rifc^iO- + am k,k+i)- in particular for a = —1 we get 

JV-l 

II( 1 -^+i) = E(- 1 ) s(J)m ^ 
fe=i j 



Now apply this with rrij^ = x(Nt (( a j 5 &&)) = 0)- These satisfy the collapsing products 
almost surely under the probability P^. The Pfafhan (|20p shows that J?£[jnj] = Pf (F| j) 
and so 



P£[iVt(4) >0for fc = l,...,2n-l] = Eg 



2n-l 



(1 - m ktk +{) 

fc=l 

= [roj] 
J 

= ^(-l)^Pf(F| J ) 
J 

= Pf(J-F). 

We may apply the same argument for the annihilating case taking m^ k = (— l) Ar *(( a i' a fe)) i 
where 1 — m^fc = 2\{N t ((aj, a k )) is odd), to find (j22j) . 

For (p3|) we have Pf(02n) = 1 and Pf(02n|j) = unless 02 n \j is a copy of Oi m for 
some m £ {0, 1, . . . , n}. This occurs either if J is empty or if J is of the form 

Ji = {2k 1 -l,2k 1 ,2k 2 -l,2k 2 ,...,2k m -l,2k m } 

for some 1 < k\ < . . . < k m < n. (36) 

Then formula (114ft specializes to 

p f (o-A) = ^(-i)i Ji i/ 2 Pf(AU) 
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where the sum is over all J\ of the form in (|36p (including the empty set). We use 
another combinatorial identity, also straightforward by induction on n: 

n 

n( 1 -^-i, 2fe ) = ^(-l) |Jl|/2 mj 1 
k=l Ji 

where the sum is over all J\ of the form in (|36p (including the empty set). Arguing as 
in the previous example leads to ([23]) . 

For PH) one has Pf(02 n ) = and Pf(02«|j) = unless 02 n \j is a copy of 2m for 
some m € {0, 1, . . . , n — 1}. This occurs either if J is empty or if J is of the form 

J 2 = {2ki,2k 1 + l,2k 2 ,2k 2 + l,...,2k m ,2k m + l} 

for some 1 < ki < . . . < k m < n — 1. (37) 

Then formula (|14[) specializes to 

Pf(,4-0) = ^(-l)^/ 2 Pf(A| c j2 ) 

where the sum is over all J 2 of the form in (j37j) (including the empty set). The required 
combinatorial identity is 

n—l ra 

IJ (1 - m 2k)2k+1 ) Yl m 2k -i, 2 k = ^(-l) |J2|/2 mj 2 
fc=i fc=i j 2 

where the sum is over all J 2 of the form in (j37[) (including the empty set). Arguing as 
in the previous examples leads to (fM|h 

4.3 Details on distributional derivatives 

The derivation of the kernel K in section 13.31 and also the Pfafhan (|27p . use formal 
differentiation that can be made precise by using distributional derivatives. Consider 
first ()27p . For fi = a locally finite point measure with disjoint atoms, one has the 

distributional derivative on V 2n 

d Xl --- d X2n x(v(xk, x k +i) = for k = 1, 2, . . . , 2n - 1) 

= {-l) n x(K x k,x k+1 ) = for k = 1,3,5. .. , 2n - l)n(dxx) .. .fi(dx 2n ). (38) 
We illustrate how to check this by showing that, in the distributional sense on V 2 , 
dxX(Kx,y) = 0) = x{^ix,y) = 0)fi(dx)dy. 
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Indeed, if / is smooth and compactly supported in V 2 then 
f{x, y)x{^{x, v) = 0)fi{dx) dy 
y2 / f{xi,y)x{^{zi,y) = 0, z, L < y) dy 

J 2 d x f{x, y) ^2 X(^(zi,y) = 0, x < Zi < y)^j dx dy 
d x f{x,y)x(v(x,y) > 0)dxdy 



d x f(x, y)x{n(x, y) = 0) dx dy 



since at most one term in the sum over i is non-zero. Iterating such calculations leads 
to ([38]) , Then for smooth / compactly supported in V 2n , 



V 2 n 



f(xi, . . .,X2n)Pt(xi, . . .,x 2n )dx 1 ...dx 



2/i 



E. 



c 



\ 2n 



f(xi, x 2n ) X (N t (I k ) = for k = 1, 3, . . . , 2n - l)N t (d Xl ) . . . N t (dx 2n ) 



n jpC 

00 



/ d Xl ... d X2n f(xi, . . . ,x 2n )x(Nt(h) = for k = 1,2, ... ,2n- l)dxi ---dx 

UV 2 n 

(-l) n [ d xl ...d X2 J(x 1 ,...,x 2n )Pi(F(t- 1 / 2 (x j -x i ))) d Xl ...dx 2n 

JV2n V ' 

(ATrt 2 )^ 2 [ f(x U X 2n )Pi {^{t- 1/2 {Xj - Xi ))) d Xl ... dx 2n . 



v 2n 

In the last step we have passed the derivatives onto the Pfaffian, which is smooth since 
F is smooth, and used F"(x) = (47r)~ 1 ^ 2 (f>(x). 

The argument for the kernel K is similar. Fix x 2 < X4 < . . . < x 2n and consider the 
open set V = {(x±, x 3 , . . . , x 2n -\) : x\ < x 3 < . . . < x 2n -\}. Then, as above, in the 
distributional sense on V 

d Xl d X3 ■ ■ ■ dx^^xi^Xk, x k +i) = for k = 1, 3, . . . , 2n - 1) 
= xiK x k,Xk+i) = for k = 1,3. . . ,2n- l)/i(dxi) fi(dx 3 ) . . . fJ,(dx 2n -i). 

Then for smooth / compactly supported in V, with Q, = {N t (Ik) = for k = 1, 3, . . . , 2n — 1}, 



d Xl d X3 . . ■d X2n _ 1 f(x 1 ,x 3 , . . .,x 2 n-i)x(fyd%i dx 3 ... dx 2n -i 



E oo 


\ f(xi,X 3 , 

Jv 


. . . , x 2 






/ 9 Xl 
Jv 




(-1)" / d xl d X3 .. 
Jv 




Pf 


[Fit'' 



Pf ( F(t~ ' (xj - x^) : 1 < i < j < 2n dxidx 3 ...dx 



t-2n-l 
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Now one can pass the derivatives onto the Pfaffian and then let x-i \, xi, X4 4, x%, . . . as 
described in section 13.31 



4.4 Details for section T3.4I 



Here we give the error estimates for the Pfaffian expansion Lemma 

The product (|30p that defines the symmetric polynomials a^ n yields a total of 2 2n 
monomials so we have the simple bound |<7fc n (y)| < 2 2n \y\ k . The expansion (|3ip must be 
iterated at most p times to derive (|32[) and this leads to to the bound 

| r 2n, P+ 2n (y) | < ( 2n 2 2n )P\y\P+ 2n - k+1 . (39) 

Using this we may bound the size of the remainder terms given in (|33|) . For example 

oo 

E i^ 1 (v)ii'«ibf , *~ l (v)i 

k,l=2n+l 

oo _. 

^ E (k _ fm - IV. I^ + '~ 2 (°)I {2n2 2n ) l+k -^- 2 \y\ l+k -^ 

k,l=2n+l ^ ''^ >' 

oo ^ 

^ ^ E 71 _ nTa - 1)! I^ + '' 2 (°)I (2n2 2 "6) /+fc - 4 "- 2 when|y|<e 

k,l=2n+l ^ '' 
oo r 

< w 2 E E ^I0 r (o)i(2n2 2 « e r 4 ™ 

r =4n \s\<r-4n 

using r = k + 1 — 2, s = k — I and > 2~ fc ~ ? 

oo 1 

< 2 4n |y| 2 E - \f(°)\ 2r(4n2 2n e) r " 4n . 

r=4n r ' 

Choosing e = e(n, 0) so that 4n2 2n e lies in the radius of convergence of <fi we obtain a 
convergent series. Similarly 



oo oo 1 

E 1^-^)1 w < E (<7 -i),a-i), ^ + '" 2 wi( 2 " 2 

i=2n+l Z=2n+1 W ; ' V ; ' 



2n^i— 2n— 1 P 



oo _. 

^ If I E ?73T ry|^- 2 (0)|(2n2 2 ™6)'- 2 "- 1 

«=2n+l ^ 
oo 

^ If I E r — T —r,^ q+l - 2 (0)\(2n2 2n ey- 2n - 1 (l + 2n) 2n 

l=2n+l ^ >' 

< C(n,<f>)\v\- 



A similar bound holds for the final term in ([33]). Combining the estimates yields the 
desired error bound on R^ 1 ^ . Moreover these bounds show the absolute convergence that 
justifies the rearrangement of the series ()29[) used in Lemma [5] provided that \y\ < e(n, </>). 
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